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ABSTRACT 

We report on the discovery of SN 2014J in the nearby galaxy M 82. Given its 
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proximity, it offers the best opportunity to date to study a thermonuclear super- 
nova over a wide range of the electromagnetic spectrum. Optical, near-IR and 
mid-IR observations on the rising lightcurve, orchestrated by the intermediate 
Palomar Transient Factory (iPTF), show that SN 2014J is a spectroscopically 
normal Type la supernova, albeit exhibiting high-velocity features in its spec- 
trum and heavily reddened by dust in the host galaxy. Our earliest detections 
start just hours after the fitted time of explosion. We use high-resolution optical 
spectroscopy to analyze the dense intervening material and do not detect any 
evolution in the resolved absorption features during the lightcurve rise. Similarly 
to other highly reddened Type la supernovae, a low value of total-to-selective 
extinction, i?y< 2, provides the best match to our observations. We also study 
pre-explosion optical and near-IR images from HST with special emphasis on the 
sources nearest to the SN location. 

Subject headings: supernovae: individual(SN 2014J) — galaxies: individual (Messier 
82) — dust, extinction 


Introduction 


Type la supernovae (SNe la) are among the most luminous transient eve nts at opt ical 
wavelengths and extremely valuable tools to measure cosmological distances, see lGoobar & Leibundgut 
12 01111 for a recent review. Yet, SNe la close enough to allow for detailed scrutiny of their 
physical properties are ve ry rare, especially i n a galaxy like M 82 , the host of several recent 
core-collapse supern ovae flMattila et al.ll2013t iGendre et al.ll2013l) . At an estimated distance 
to M 82 of 3.5 Mpc I Dalcanton et al. 2009 1. SN 2014J is th e closest identified Type la S N in 
several decades, possibly r ivaled by SN 1972E in NGC 5253 f Ardeberg &; de Groot 1973 1 and 
SN 1986G in NGC 5128 ( Phillips et al. 1987 1. Thus, SN 2014J is exceptionally well-suited 
for follow-up observations in a wide range of wavelengths, from radio to gamma-rays. These 
have the potential to yield transformational new clues into the progenitor systems of SNe la, 
as well as the detailed properties of dust along the line of sight, key astrophysical unknowns 
for the study of the accelerated expansion of the universe. 

There is strong evidence that SNe la ari se from thermonuclear explosions of carbon- 
oxygen white dwarfs (WD) in binary systems (INugent et al.ll201ll: iBloom et al.ll2012l ) . How- 
ever, the nature of the second star remains unclear. For a long time, the preferred scenario 
was the single degenerate (SD) model fl Whelan & Iben 1973), where a WD accrets mass 
from a hydrogen or helium rich donor star, thus becoming unstable while approaching the 
Chandrasekhar mass. The double-degenerate (DD) model involving the merger of two WDs 
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(Tutukov & Yungelson 

1981: 

Iben & Tutukov 

1984; 

Webbink 

1984 

) has gained 

)• 

observational support in recent years, see e.g., 

Wang & Han 

(2012 


considerable 


In this work, we search for potential signatures of a SD progenitor system, such as 
variable Na D lines, precursor nova eruptions, features in the early lightcurve, radio emission, 
or a coincident source in pre-explosion in HST images. 


2. Discovery and classification 


SN 2014J was discovered by Fossey et ah (2014) in BVR - band images of M 82 obtained 
on 2014 January 21.81 UT at UCL’s University of London Observatory. We have performed 
image subtractions using pre-explosion data from the Palomar P60 telescope as reference, 
calibrated with nearby stars listed in the APASS catalog] yielding a discovery magnitude of 
R = 10.99 ± 0.03 mag. The discovery image (S. Fossey, private communication) and the P60 
reference image, as well as the difference between the two are shown in Fig. [0 along with 
the pre-explosion HST images (GO: 11360, PI: R. O’Conell; GO: 10776, PI: M. Mountain). 
The relative position of SN 2014J with respect to neighboring stars (middle panel) was 
established using multiple sh ort expos ures in Ji-band with adaptive optics and the NIRC2 
wide camera at Keck ( Tendulkar et al. 2014). In Section [5] we present a detailed analysis of 
the pre-explosion data. 

A classification spectrum was obtained by the iPTF team on January 22.30 with the 
Dual Imaging Spectrograph on the ARC 3.5m telescope ([.Cap et ah 2014), and in the near-IR 
using the MOSFIRE instrument at Keck. The combined spectra are shown in the bottom 
panel of Fig. Q3 while the photometry collected so far is displayed in Fig. [21 The object 
shows characteristic spectral features associated with SNe la, e.g., similar to SN 201 lfe 
(IPereira et al.l 120131: iHsiao et al.l 120131) . However, the steep attenuation of the spectrum at 
short wavelengths is indicative of unusually large extinction by dust in the line of sight. A 
good match to the overall SED is found invoking a pronounced color excess, E(B — U)host ~ 
1.2 mag, in addition to Galactic reddening, E(B — U) MW = 0.14 mag (jSchlaflv & Finkbeiner 
201l|), as shown in Fig. [3l For the comparison, the spectra o f SN 2011fe were artificially 
reddened assuming a Milky- Way type extinction law flCardelli et al.l 19891 1. where both the 
color excess and Ry were allowed to vary freely. The spectrum favors a low value of the 
total-to-sele c tive e xtinction, R v <2, as also suggested by spectropolarimetry observations by 
Pat at et al.l (2014). Low values of R y are not unusual in SNe la, especially in the cases of 
high extinction, see e.g. Nobili & Goobar f2008). 


1 http://www. aavso.org/apass 
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3. The iPTF-led multi-wavelength monitoring of SN 2014J in M 82 


As a part of its continuous survey of the sky in the search for transients, iPTF has 
monitored M 82 since October 2009, with nearly daily cadence over the several months 
each year when M 82 is visible from Palomar. The most recent campaign started on 2013, 
November 28. For the periods around full moon (e.g., around the time SN 2014J exploded), 
not well suited for transient searches, Ha narrow-band imaging was conducted. 


The current best fit of the time of explosion, to, was reported by KAIT flZheng et al. 


2 01 4) to be January 14.72 UT (± 0.2 days). Upon later scrutiny of the pre-discovery P48 
data, the supernova was found in several observations from the iPTF Ha narrow-band 
survey, starting just hours after the fitted to- We find a relative flux in crease from Jan uar y 
15.18 to January 16.18 of 1.6 mag, consistent with the “method 2 fit” in lZheng et ah, I (2 0141 ) . 


The SN is also prominent on R-band photometry from the P48 prior to January 21 shown 
in Fig. [21 but remained undetected by our automated software due to pixel saturation. 

Through an iPTF-led effort, involving _also the Las Cumbres Observatory Global Tele- 


scope (LCOGT) n etwork flBrown et al. 


Abu Observatory (iVenkat a raman et al. 


2013), the Nordic Optical Telescope and the Mount 


2 0141 ). we have secured optical, near-IR and mid-IR 


lightcurves carefully monitoring the rise of the supernova, as shown in Fig. [21 The 4.5/irn 
observations were taken under the Spitzer InfraRed Intensive Transients Survey (SPIRITS; 
PI: Kasliwal). 

The spectra shown in Fig. [3] are consistent with those from a normal Type l a expl osion, 
similar to e.g., SN 2011fe, but reddened following a CCM-law flCardelli et al.l 1989]) with 
E(B — V) ~ 1.2 mag and Ry = 1.3 — 2, in addition to G alactic reddening. Fig. [21 also 
shows lightcurve fits using the SNooPy fitter ([Burns et al. 201111 of the photometric data 
prior to maximum brightness. Best fits are found for E(B — U)host = 1-22 ± 0.05 mag 
and Ry = 1.4 ± 0.15. We expect the accuracy of the fitted parameters to improve, as the 
lightcurve shape estimate will profit from the decreasing part of the SN lightcurve. However, 
the available data clearly puts Ry well below the Galactic average value, Ry = 3.1. 

We also obtained two 1800 s high-resolution (R = 40000) spectra with SOPHIE at 
Observatorie Haute-Provence on January 26.0 and January 28.0. Further, two 1800 s spectra 
(R = 67000) were obtained with the Fibre-fed Echcllc Spectrograph (FIES) on January 27.3 
and another on February 1.0 with the Nordic Optical Telescope. 

All spectra reveal deep multiple component Na I D absorption and diffuse int erst ellar 
bands (DIBs), including AA 5780,5797,6284, and 6614, also reported by Cox et al. ( 2014 1 and 
Kotakl (2Q14). The SOPHIE spectra further contain well resolved Ca II H & K with features 
matching those of the Na I D lines, shown in the co-added spectrum in FigQ] We have not 
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detected any significant time evolution for any of the resolved components of the Na I D 
doublet over the four epochs (at the ~10% level for 3cr), thus motivating the combination 
of the spectra. 


Following the procedure outlined by Phillips jet ah (120131 ). we measure the equivalent 
width (EW) of the A5780 DIB to derive an independent estimate of host galaxy extinction 
for SN 2014J. We find EW(5780) = 0.48 ±0.01 A corresponding to Ay ost = 2.5 ± 1.3 mag. 


Given the low recession velocity of M 82, it is difficult to separate the contribution from 
the Milky- Way and the SN host galaxy absorption. However, the availabilit y of H I data from 
the LAB survejl^] in the direction of M 82 (see inset panel in Fig. [4] Kalberla et al. (2005)) 
clearly indicates which features are Galactic. Hence, all the absorption features redshifted 
with respect to the Milky Way are due to intervening material in M 82. 


4. Spectral modeling 


In Fig. |3]we present a time-series spectral comparison between SN 2011fe and SN 2014J 
starting roughly 12 days before peak brightness. The SNe are remarkably similar in their 
spectral evolution. The main differences seen are that the overall velocities are higher in 
SN 2014J (see the inset Si II velocity plot) and there is a strong signature of high-velocity 
Si II and Ca II in this supernova. 


To further investigate these differences we carried out a set of SYNAPPS (Tho ma s et al. 


2011) fits to these two supern ovae as we ll as to SN 2005cf, which was distinct in its pervasive 


high-velocity features (iWang et all 20091 ). We present the results in Fig. 0 In our fits to 
the red-side of the optical spectra we have employed the ions: C II, O I, Mg II, Si II and 
Ca II with the latter two having both photospheric and high-velocity components. We see 
that SN 2014J more closely resembles SN 2005cf with respect to the high-velocity Si II and 
Ca II features which extend easily over the range of 20,000-30,000 kms -1 . Unlike either 
SN 2005cf or SN 2011fe, C II and O I are absent at this phase in SN 2014J. We searched 
for the presence of C I A 1.0693 /jrn line in our NIR spectrum, but the signal-to-noi se is too 
low to do a meaningful fit comparable to that done for SN 2011fe (IHsiao et ah 20131) . We do 
note that due to the stronger than average Mg II features seen in the optical for SN 2014J, 
this analysis may be more challenging, even with higher quality spectra, as the presence of 
Mg II A 1.0927 jum should be quite strong as well. 


Several explanations for the origin of the high-velocity features have been presented, 


2 http://www. astro, uni-bonn.de/en/download/data/lab-survey/ 
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from density enhancements via swept up CSM ( 

Gerardv et al. 

2004: 

Tanaka et al. 

20061) to 

mixing or more complete burning in the outer layers of the supernova i 

Mazzali et al. 

2005aJ|b|) 


What is clear is that the features 
do offer a unique diagnostic for understanding properties of the progenitor system and/or 
the explosion mechanism and correlations betw een the strength of these features and th e 
underlying colors and lightcurves of the SNe la flChildress et al.l I2014J : iMaeuire et al.ll2013h . 


SN 2014J is among a class of SNe la where high-velocity features are present yet little 
to no evidence for C If ex ists eve n in very early spectra (see the broad line or high-velocity 
gradient examples in Parrent et al. 12 0111 )1. Since extensive UV data from HST exists for 
both SNe 2005cf and 201 lfe, it will be interesting to see which of these supernovae SN 2014J 
most closely matches with respect to both the color and luminosity evolution. 


5. The quest for the progenitor system 


M 82 has been extensively imaged by HST, thus it is possible to study the environment 
of the SN prior to the explosion. Because of the large attenuation due to dust in the line of 
sight, we concentrate on the NIR bands. We perform aperture photometry on the nearest 
sources to SN 2014J, shown in Fig.[D The closest object (yellow circle in Fig. [I]) falls 0.2" from 
the current best estimate of the SN location (RA = 9 / i 55 Tn 42 s .217(l), Dec = 69°40'26".56(4) 
in J2000 coordinates with respect to the HST image, Tendulkar et ah 2014), corresponding 
to a 4cr spatial offset. For this source, we measure AB magnitudes of FTlObF = 21.4 ± 0.4, 
F128N = 21.9±0.4, JT60JF = 20.6±0.4, F1641V = 21.2±0.4. The error is dominated by the 
uncertain background subtraction as a result of source confusion. The F110 — F160W color 
is typical of other sources near this position. At a distance of 3.5 Mpc, the corresponding 
absolute magnitude of the nearest resolved object is Jab Hab ~ -7 mag (A H < Aj < 0.4 
mag). This source could represent a stellar cluster, a grouping of unrelated objects or a 
region of relatively low dust attenuation. 

Next, we consider the possibility that the source is a donor star in the SD scenario. 
The derived luminosity would then suggest a very luminous red supergiant. However, in 
the case where a SN originates from a system with CSM created by a mass-loss from a 
donor star, the int erac tion b etwe en the SN ejecta and the CSM is expected to give rise 
to radio emission 


Chevalier 


19821 1988th The r adio null-detections on January 23 and 


24 (jChomiuk et al.ll2014j; IChandler &; Marvill 120141) can therefore be used to derive an upper 
limit on the mass-loss rate. Adopting a SN shock-wave velocity of 3 x^lO 4 km s^ 1 (ab out tw ice 


the Si 11 velocity, see Fig. [3]), and similar parameters as assumed by Horesh et ah (2012) for 
SN2011fe, the upper limit on the mass-loss rate is M < 7 x 10 _9 (tc/100kms _1 ) M Q yr _1 , 


where w is the mass- loss wind vel o city. The upper limit is comparable to the ones obtained 
for SN 2011fe (Horesh ct al. 2012 ; Chomiuk et ah 2012 b Given these tight limits and the 
spatial displacement, we conclude that the closest resolved source in the pre-explosion HST 
images is unlikely to be a donor star. 


Finally, we note that the FTIOW — FT60JF color map shown in the middle right panel 
of Fig [□ suggests that SN 2014J is at the edge of a dust patch, about 4 pc in projected size. 
Light echoes may thus be expected for this supernova. 


We have searched for possible nova outbursts in the historic P48 i?-band data covering 
a period of about 1500 days prior to the detection of SN 2014J. By binning the data in 
bins of 15 daya, we dcy not find any excess larger than 4cr (calculated using the bootstrap 
technique, see Efron (1982))- Our limiting magnitude is R > 19.5 mag for a total time span 
of 1000 days, and R > 20.25 mag for more than 765 days in this 1500 days time window. 
Assuming Ar = 2 mag, compatible with the extinction we estimate based on the supernova 
colors, this corresponds to absolute magnitudes Mr = —10.2 and —9.5 respectively at the 
dist ance o f M 82. However, given the uncertainties on the properties of recurrent novae, see 
e.g. ITang et all (20141). we refrain from drawing firm conclusions against the possibility of 
recurrent novae preceding SN 2014J based on these non-detections. 


6. Summary and conclusions 

The discovery of SN 2014J presents us with a unique opportunity to explore the physics 
of Type la SNe and the line-of-sight effects due to intervening matter. Further understanding 
in these areas is of utmost importance for the use of SNe la in cosmology. The early data 
from P48, starting as early as only hours from the explosion, and the multi-wavelength 
follow-up by the iPTF team covers an important range of the available windows in the 
electromagnetic spectrum. Just as the lightcurve reaches its maximum, we have learned 
that the SN has suffered non-standard extinction. We have searched for, but not detected, 
any time variation in our high-resolution spectra of the Na I D doublet. Similarly, we do 
not detect any pre-explosion activity in the ~1500 days of P48 monitoring. In a study of 
pre-explosion HST images in the near-IR, the nearest resolved source is found 0.2" away 
from the SN location. The source brightness and offset from the SN makes it unlikely as a 
donor star in a single-degenerate scenario. 

Further, we make a first study of the spectral features of SN 2014J and find that it 
exhibits high-velocity features from intermediate mass material but lacks C and O often 
seen in very early spectra. Otherwise, it is a very similar to several well-studied normal SNe 
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Discovery image, Jan 21 .81 


Reference image 


Subtraction 



A 

Keck NIRC2 AO, Jan 22.38 



-I n 


HST F8 1 4W/F 1 1 0W/F 1 60W 


SN 201 4J, ARC DIS +Keck MOSFIRE, Jan 22 - 
De-reddened with E(B-V) host = 1.2 ( R v = 1.7) 
and E(B-V) msn = 0.14 {R v = 3.1) 

SN 201 Ife - 



5000 


10000 

Restframe wavelength, A rest [A] 
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Fig. 1. — Top panels: 5' x 5' sections of the discovery ( Fossev et al.l 2014). reference and the 
subtraction images. Middle panels (from left to right): 5" x 5" section of the Keck NIRC2 
adaptive optics image used to match the SN coordinates (red circle) to the surrounding 
sources (black circles). The nearest resolved object (yellow circle) in the pre-explosion HST 
composite images (F815fF/F110fF/F160kF) is offset by 0.2" from the best estimate of the 
SN position. The middle-right panel shows a color map (FT10PF — FT60PF) indicating the 
large scale structures, probably due to patches of dust. The bottom panel shows the first 
optical and NIR spectrum of SN 201 4J and a comparis on to a combined spectrum from 
SN 2011fe by Pereira et al. ( 2013 ) and Hsiao et ah ; ( 2013 ). described in the text. 
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Fig. 2. — Lightcurves showing the rise of SN 2014J, along with SNooPy fits described in 
the text. The first P48 H® 56 and _H^ n arrow bands detections are shown (black circles 
and inset image), S-corrected (IStr itzi nger et ah 2002) to /I'-band. Due to lack of accurate 
absolute calibration for the H a filters, a common offse t was ap plied to co nnect with the 
data-points in the pre-discovery light curve presented in iZheng et all (2014). We also show 
their two fits of t 0 , suggesting our first detection could be within ~5 hours from the onset 
of the supernova. 
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Fig. 3. — Pre-max spectroscopic follow-up of SN 2014J, starting 7.6 days from estimated 
supernova onset using the ARC 3.5m DIS, LCOGT FLOYDS, NOT ALFOSC and P200 
DBSP spectrographs. Spectra of SN 2011fe ( Pereira et ah 20131 ) at similar epochs (gray 
lines, reddened by E(B — R)host = 1-2 mag with R v = 1.7 and E(B — R)mw = 0.14 mag 
with Ry = 3.1 are shown for comparison. The inset panel shows that SN 20 MJ has higher 
Sill velocity than SN 201 lfe and a steeper velocity gradient (iBenetti et ah 2005). 
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Fig. 4. — The normalized Na I D doublet is plotted with solid vertical lines indicating the 
Galactic rest, fra m e wavelengths and the dashed vertical lines correspond to the mean velocity 
of M 82 (203 kms -1 ). The inset panel shows the velocity distribution of Ca II ff & K (green 
and red lines) and the A5780 DIB (black). Additionally,_the_H I A21_cm emission spectrum 
of the line of sight of M 82 from the LAB survey ( Kalberla et ah 2005) is plotted (dashed 
gray line). The features of the Na I D and Ca II H & K, at approximately -50 and 0 kms -1 
with respect to the Galactic rest frame can be attributed to absorption in the Milky Way. 
The main Na I D absorption features, saturated between 74-135 kms -1 , originate in M 82. 
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Fig. 5. — Early spectra of SNe 2011fe, 2005cf and 2014J, concentrating on the red-side of 
the optical spectra, along with SYNAPSS ( Thomas et al. 2011) fits to each (dashed lines). 
In these fits we have employed the ions: C II, O I, Mg II, Si II and Ca II with the latter 
two having both photospheric and high-velocity components. These fits highlight the major 
difference between these supernovae. SN 2014J shows no sign of either C II or O I in this 
part of the spectrum and, unlike 201 lfe, quite strong high-velocity components of Si II, and 
Ca II extending well beyond 20,000 kms -1 . In addition, the Mg II photospheric feature in 
SN 2014J is stronger as well at this phase. SN 2005cf shows similar high-velocity Si II and 
Ca II to SN 2014J, but differs from this supernova as photospheric C II is clearly seen at this 
phase and our SYNAPPS fits also favor O I at this phase. This suggests that the nuclear 
burning in the outer layers of SN 2014J was more complete than that of SN 20 life and the 
lack of C II compared to SN 2005cf may impl y the same or perhaps a viewing angle effect 
due to an off-center detonation (IParrent et ah 2011). 


